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THERMAL  EVAPORATION  IN  BLAST  WAVES:  A  MODEL  FOR 
ANEURISM  FORMATION  IN  THE  NRL  LASER-HANE 
SIMULATION  EXPERIMENT 


I.  INTRODUCTION 

During  the  past  two  years  the  Plasma  Physics  Division  at  the  Naval 
Research  Laboratory  has  undertaken  a  program  to  develop  a  set  of 

experiments  which  would  simulate  certain  physical  processes  that  are 
characteristic  of  a  high  altitude  nuclear  event  (HANE).  The  program  is 
under  the  auspices  of  DNA's  Division  of  Atmospheric  Effects  and  has  the 

primary  goal  of  a  better  understanding  of  the  degradation  of  radar  and 
communication  transmissions  through  an  atmosphere  seeded  by  nuclear 

bursts.  At  the  present  time  the  experiment  is  initiated  by  irradiating  a 
thin  target  foil  with  a  neodymium  laser.  The  subsequent  rapid  heating  and 
disintegration  of  the  target  is  followed  by  an  expanding  shock  wave  once 
the  target  debris  couples  to  the  background  gas.  Typical  experimental 

parameters  used  thus  far  consist  of  an  Al-foil  less  than  1  mm  in  diameter 
with  a  total  mass  of  several  tenths  of  a  Ugm,  a  3.4  nsec  laser  pulse  of 
total  energy  4-400  Joules,  and  a  mostly  background  gas  of  pressure 

0. 001-10. 0  Torr  threaded  by  a  uniform  magnetic  field  of  either  0  or  800 
Gauss  strength. 1 

An  important  diagnostic  has  been  dark-field  shadowgraphy  which  gives  a 

visual  impression  of  the  coupling  region  (or  blast  wave)  as  it  expands. 

This  process  picks  out  steep  density  gradients  in  the  image  plane  and 

2 

permits  multiple-time  recordings  on  a  single  photograph.  Stamper  et  al. 
present  several  photographs  which  depict  the  coupling  region  at  52  and  96 
nsec  after  the  initial  laser  pulse,  when  the  blast  wave  has  expanded  to 
about  1  cm  from  the  target  position.  An  unusual  feature  seen  in  the  high 
ambient  pressure  shots  is  the  development  of  an  aneurism  in  the  coupling 
shell:  an  outward  protrusion  from  the  nearly  spherical  blast  wave  toward 

the  direction  of  the  laser.  The  purpose  of  this  report  is  to  present  a 
physical  mechanism  which  may  account  for  the  development  of  this 

aneurism.  Since  this  feature  seems  to  ocur  only  in  shots  with  a  high 
ambient  pressure  (several  Torr),  the  present  theoretical  model  may  have 
relevance  for  HANE  phenomena  at  altitudes  of  100  to  200  km. 

Manuscript  approved  July  11,  1984. 
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II.  OVERVIEW  OF  THE  PHYSICAL  MECHANISM 


In  this  section  we  develop  a  schematic  picture  of  the  proposed 
physical  mechanism  responsible  for  the  aneurism.  A  detailed  account  of  the 
specific  model  equations  will  be  given  in  the  following  section.  Ripin  et 
al .  have  used  the  dark-field  shadowgraphs  to  show  that  the  radius  of  the 
coupling  shell,  R,  expands  as  an  adiabatic  blast  wave.  The  Taylor-von 
Neumann-Sedov  self-similar  solution  shows  that  R  ■  £  (E  /p  where 
Eq  is  the  total  laser  energy,  pQ  is  the  ambient  density,  t  is  the  time 
since  the  laser  pulse  and  is  a  constant  of  order  unity.  The  basic 
physics  of  the  blast  wave  can  be  readily  found  by  assuming  that  all  of  the 
swept-up  mass  resides  in  a  thin  shell  moving  with  the  leading  strong  shock 
front  while  being  pushed  by  the  pressure  inside  the  enclosed,  hot  cavity^. 

Basically,  the  proposed  model  complicates  this  simple  picture  by 
including  a  thermal  candution  front  which  forms  on  the  inside  edge  of  a 
small  segment  of  the  swept-up  shell.  The  remainder  of  the  shell  is 
thermally  insulated  due  to  the  magnetic  fields  which  are  generated  at  the 
target  during  laser  burning  and  are  subsequently  carried  out  with  the 
expanding  plasma.  Within  the  conduction  front  the  electron  heat  flux  from 
the  hot  cavity  to  the  cooler  shell  is  balanced  by  mass  evaporation  off  the 
shell  with  the  flow  directed  into  the  cavity.  The  mass  unloading  of  the 
evaporating  segment  results  in  an  acceleration  relative  to  the  insulated 
segment,  assuming  the  cavity  pressure  Is  nearly  uniform.  Furthermore,  the 
evaporating  segment  is  accelerated  by  a  rocket-like  effect;  the  inward 
momentum  carried  off  by  the  evaporative  mass  flow  is  balanced  by  forwad 
moving  pressure  waves  driven  into  the  swept-up  shell.  Of  course,  as  the 
evaporating  segment  is  accelerated  forward,  it  sweeps  up  relatively  more 
material  and  suffers  an  Increase  in  the  opposing  ram  pressure.  Under  the 
proper  conditions  we  will  show  that  this  evaporating  segment  moves  ahead  of 
the  remaining  shell  resulting  in  an  aneurism-like  feature. 

A.  Thermal  Conduction 

We  first  justify  the  model  by  estimating  some  relevant  physical 
quantities.  The  effect  of  thermal  conduction  is  central  in  the  proposed 
model,  but  how  Important  is  thermal  conduction  in  the  experiment?  Consider 
the  total  energy  equation  for  the  cavity, 
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ev”^/2.  xhe  constant  K  was  calculted  from  Braginskii^  for  a  Coulomb 
logarithm  of  10  and  a  charge  per  ion  in  the  cavity  Zc  ■  5.  Let  L^.  be  the 
temperature  scale  height,  L  *  T/|$T|,  then 
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where  Tc  is  a  characteristic  cavity  temperature.  The  300  eV  estimate  for 

Tc  was  based  on  the  following  considerations.  The  high  sound  speed  within 

the  cavity  implies  a  nearly  uniform  cavity  pressure  pc  and  the  conservation 

of  total  energy,  together  with  an  estimate  of  the  postshock  gas  velocity, 

leads  to  the  relation  p  -  p  /2  ,  where  p  is  the  shell  pressure^.  Hence, 
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where  Zg  is  the  mean  charge  per  ion  in  the  shell.  The  mean  mass  per 

particle,  u,  has  been  replaced  by  m^/(  1  +  Z),  with  the  mass  of  a 

Nitrogen  ion.  The  measured*5  values  at  the  time  of  aneurism  formation 

are  T  **  15eV,  Z  '  3  and  Y  ■*  1.2  implying  p  /p  '  11,  where  P  is  the 

s  s  s  so  o 

ambient  density.  Since  most  of  the  swept  up  mass  resides  in  the 

shell,  P  /P  »  11  .  In  our  calculations  we  will  consider  the  set  of 

s  c 

values  T  /T  *  (10,20,30)  with  the  middle  one  corresponding  T  *  300  eV 
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Let  us  now  estimate  the  enthalpy  flux  Yvp/(Y  -  1)  in  the  cavity.  The 
3 

kinetic  energy  flux  pv  /2  in  the  energy  equation  (1)  is  negligible  due  to 

the  low  cavity  density.  Again  using  p  -  p  /2 
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where  y  =  y  =  5/3  and  3  x  10  cm  and  10  cm  sec  are  a  typical  electron 

c  6 
density  and  shell  velocity,  respectively,  for  the  high  pressure  shots  . 

Clearly  the  approximately  measured  enthalpy  flux,  Eq.  (4),  is  much  less 

than  the  estimated  heat  flux,  Eq.  (2)  for  any  reasonable  L^,  (<lcm)  .  A 

similar  result  holds,  though  to  a  lesser  degree,  even  if  the  heat  flux  is 
saturated.  Saturation  of  the  heat  flux  will  be  discussed  more  fully  in 
section  III.  Thus  the  expansion  of  the  cavity  should  be  more  nearly 

isothermal  rather  than  adiabatic  as  in  the  Taylor-von  Neumann-Sedov  self¬ 
similar  solution. 

3.  Magnetic  Fields 

Since  the  thermal  conductivity  perpendicular  to  the  magnetic  field  can 
be  smaller  than  along  the  field  it  is  important  to  Investigate  the  topology 
and  strength  of  the  magnetic  field  in  the  experiment.  Consider  the  case 

where  there  is  no  ambient  magnetic  field.  According  to  Kacenjar  et  al.^  a 
magnetic  field  of  100  to  200  Gauss  was  still  detected  as  the  shell  crossed 
the  first  magnetic  probe  located  1  cm  from  the  symmetry  axis  and  R  *  1  cm 
from  the  target.  A  second  measurement  at  the  same  offset,  but  R  *  3  cm 
from  the  target  found  that  the  field  strength  near  the  shell  had  decreased 
to  less  than  20  Gauss.  The  orientation  of  the  magnetic  probe  indicated 
that  the  field  was  azimuthal,  i.e.,  field  lines  circled  the  symmetry 
axis.  The  presence  of  this  field  agrees  with  the  theoretical  model  of 
magnetic  field  generation  through  the  thermoelectric  effect  at  the  target 
surface  during  laser  burning®"'*®.  If  this  self-generated  magnetic  field  is 
carried  outward  by  the  debris,  it  would  decrease  as  1/R  ,  consistent  with 
the  above  measurements.  This  suggests  that  the  magnetic  field  strength 
is  >  1  kilogauss  when  R  >  1  cm,  where  the  aneurism  appears  to  start 
growing.  Actually,  this  estimate  is  a  lower  limit  since  the  magnetic  probe 
may  not  resolve  the  magnetic  field  structure  and  the  probe,  which  is 
encased  in  quartz  tube  to  prevent  space  charging  effects,  requires  the 
diffusion  of  the  field  out  of  the  gas  for  a  positive  detection. 

Thus,  we  adopt  the  estimate  of  1  kG  for  the  off-axis  self  generated 
magnetic  field  strength  when  R  <  1  cm  and  an  azimuthal  orientation.  The 
topology  of  such  a  field  requires  that  it  vanish  along  the  symmetry  axis. 
(A  schematic  diagram  of  the  evolution  of  the  blast  wave  is  shown  in  Figure 
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1).  In  the  off-axis  portion  of  the  shell  the  heat  flux  is  perpendicular  to 
the  field  and  is  thereby  reduced  relative  to  that  in  the  segment  of  the 


shell  along  the 

symmetry  axis  where 

there 

is 

no  field. 

We  can 

quantitatively 

estimate  this 

effect 

by  using  the 

approximation  ■* 

<„/  [1  +  (u>t)  2  ] ,  where 
s  e 

to  is 
e 

the 

electron 

cyclotron 

frequency  and  is  the  electron  collision  time.  The  correct  relation  is 
more  complicated^,  but  the  above  simple  one  does  have  the  proper  dependence 
in  the  extreme  limits.  Using  the  previous  estimates  for  the  shell 
parameters  one  finds 
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which  means  that  in  the  shell  <,  <,  . 

1  9 

However,  the  field  cannot  remain  in  the  shell.  The  magetic  Reynolds 
number  of  the  shell  is 
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where  <J^  is  the  perpendicular  electrical  conductivity  and  is  the 

magnetic  length  scale.  Using  the  shell  thickness  Lg  <  0.03  cm 

for  Lv,  and  Z  -  3,  one  finds  R„  <  1  .  This  implies  that  the  field  will 

Ms  P1,S 

diffuse  out  of  the  shell.  Ahead  of  the  shell  R^  remains  <  1  due  to  the 
small  temperature  and  velocity,  and  a  broad  zone  with  a  gentle  compression 
of  the  magnetic  field  is  formed.  On  the  other  hand,  as  the  field  diffuses 
back  into  the  cavity  it  experiences  a  rapid  increase  in  the  temperature. 
The  scale  length  of  the  diffusion  region  in  the  cavity  where  R^  c  <  1  is 
given  by 


t  ✓  ,  c  ir,-4  _  , 107cmsec 
L  <  1,5  x  10  Z  ( - 

M  C  V 


^  ,300  eV, 
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which  is  much  less  than  the  radius  of  the  shell.  Hence  the  field  will  stop 
diffusing  once  it  moves  back  into  the  cavity  while  it  will  readily  move 
ahead  of  the  blast  wave  (see  Figure  Ic).  This  picture  is  verified  by  the 
temporal  behavior  of  the  field  strength  as  measured  by  a  stationary  probe 
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in  the  experiments  with  an  ambient  magnetic  field  . 
the  cavity  one  finds 
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Thus  away  from  the  symmetry  axis,  the  field  resides  just  inside  the 
shell  and  <  "  ,  while  along  the  axis  B  *•  0  and  <  "  Xj  .  The 

effects  of  an  ancient  magnetic  field  will  be  briefly  discussed  in  the  last 
section. 


C.  Evaporative  Mass  Loss 


Given  that  thermal  evaporation  is  confined  to  a  small  segment  of  the 
shell  the  final  question  is  how  significant  is  the  evapaorative  mass 
loss.  We  model  the  evaporative  process  as  a  quasi-steady  thin,  i.e., 

A 

planar,  conduction  front  at  the  rear  of  the  dense  swept-up  shell.  Let  v  be 
the  gas  velocity  relative  to  the  conduction  front.  For  subsonic  flow  one 
obtains  from  Eq.  (1), 
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where  uc  »  m^/(l  +  Z^)  is  the  mean  mass  per  particle  in  the  cavity. 
Following  the  nomenclature  prior  to  Eq.  (2)  we  compare 
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where  R  is  the  velocity  of  the  shell.  This  result  means  that  evaporation 
can  significantly  reduce  the  mass  in  the  conductive  segment  of  the  shell 
resulting  in  a  protuberance  (see  Figure  Id).  A  more  realistic  estimate  of 
the  ratio  will  depend  on  the  additional  effects  of  a  staurated  heat  flux, 
an  Increase  in  the  mass  gain  rate  as  the  shell  segment  accelerates  forward, 
and  a  temporal  decrease  in  the  cavity's  temperature. 


III.  MODEL  EQUATIONS 
A.  Geometry  and  Nomenclature 

Now  that  the  general  scenario  of  the  proposed  physical  mechanism  for 

the  development  of  the  aneurism  has  been  presented,  we  will  develop  a  set 

of  more  detailed  equations  to  model  the  evaporative  effects  only.  The 

topology  and  strength  of  the  magnetic  field  is  not  calculated  here;  we 

assume  it  has  the  topology  described  in  section  II.  The  specific  geometry 

is  given  in  Figure  2.  Thermal  evaporation  occurs  within  a  sector  29  , 

while  outside  of  this  sector  the  shell  is  assumed  to  be  insulated  by  the 

azimuthal  self-generated  magnetic  field.  We  will  use  the  "thin-shell" 

approximation  with  R'  being  the  radius  of  the  insulated  shell  segment,  and 

Rg  the  radius  of  the  evaporating  segment.  Primed  quantities  will  aways 

refer  to  the  former  segment  and  a  subscript  s  denotes  shell  referenced 

values.  We  will  neglect  the  detailed  shape  of  the  protrusion  and  consider 

the  evaporating  segment  to  be  part  of  spherical  surface.  This  of  course 

greatly  simplifies  the  calculation  by  reducing  the  model  equations  to  one 

spatial  dimension;  however,  it  allows  one  to  study  the  nonlinear  physical 

processes  acting  to  generate  the  protrusion.  M  and  M  are  the  masses  of 

s  ..  s 

the  respective  sections  of  the  shell,  and  vg  and  vg  the  gas  velocities  in 
the  lab  frame.  The  shock  front  on  the  outside  edge  of  the  swept-up  shell 
will  be  taken  to  be  a  strong  adiabatic  shock.  Hence  we  can  use  the 
relations 
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s  Y  +  To  s’  »  (Y  +  1)*  s 
s  s  s 


where  Rg  *  dRg/dt,  Yg  and  are  the  specific  heat  ratio  and  mean  mass  per 
particle  in  the  shell,  respectively,  and  PQ  is  the  ambient  density. 


B.  Dimensional  Form 


In  our  model  the  total  energy  EQ  is  conserved,  i.e.,  there  are  no 
radiative  losses,  and  the  cavity  temperature  and  shell  pressure  are 
approximated  as  uniform.  The  total  energy  is  the  sum  of  the  cavity's 
internal  energy  and  the  kinetic  as  well  as  internal  energy  in  the  two  shell 
segments: 
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tt  [R'  cos“(y0)  +  Rg  sin  (^e)]pc 


+  in  v  2  +  hi'v'2  +  — i— r-4irR  2ARp  +  Ur4itR"24R'p'.(12) 

2  s  s  2  s  s  y  -1  s  *s  y  “1  s  s 
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The  kinetic  energy  of  the  cavity  material  can  be  neglected  due  to  its  low 

density  and  small  velocity.  The  presence  of  a  conduction  front  does  not 

alter  the  total  energy  equation  since  it  only  acts  to  exchange  energy  from 

internal  to  kinetic.  The  internal  energy  of  the  shell  is  typically 

negelected^,  though  actually  it  equals  the  shells  kinetc  energy.  The 

equality  follows  from  the  thin-shell  approximation  which 
2 

gives  M  »  4xR  AR  p  ,  and  Eq.  (11).  The  resulting  total  energy  equation 
s  s  s  s 

can  then  be  written  as 


2  2 
+  M  v  +  H'v' 
s  s  s  s 


The  total  momentum  equation  represents  a  balance  between  the  rate  of 
change  of  the  momentum  of  the  shell  and  the  pressure  force  supplied  by  the 
cavity.  Since  there  are  two  shell  segments  one  has 


at<v;>  ■  *"K'2c°*24«>i’c  • 


Jt(v,>  -  4-R32sln2(-j)  (Pc  +  Pevp)  . 


The  additional  pressure  Pevp  in  Eq.  (15)  is  due  to  the  rocket-like  effect 
as  mass  evaporates  off  the  shell  segment  and  flows  backward  into  the 
cavity.  Its  value  will  be  developed  when  we  discuss  the  equations 
describing  the  conduction  front. 

The  total  masses  in  the  different  segments  of  the  swept-up  shell  are 
Ms  “  T1  Rs3poco92(T9)’  (16) 
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3  s  o  2  evp 


where  is  the  total  mass  which  has  been  evaporated  from  the  conducting 

segment  of  the  shell.  Again  its  value  will  be  discussed  below. 

Let  us  pause  in  the  development  of  the  model  equations  and  consider 

the  solutions  in  the  purely  adiabatic  case,  i.e.,  0*0.  The  results 

will  be  required  as  initial  conditions  prior  to  the  onset  of  thermal 

evaporation  over  a  segment  of  the  shell.  Combining  Eqs.  (11),  (13),  (14), 

and  (16),  one  readily  finds 

R  -  C  (Vv/5  ,  ( 18) 

s  op 


5  r  75  (Ws  *  L) 

o  =  16*  [Y  +1  +  4  (y  -1)]  ’ 
s  c 


and 


_ 1  •  2 

Y  +1  Po  Rs 
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1 

7?s 


(20) 


These  results  agree  with  Ripin  et  al^.  If  the  shell's  internal  energy  had 

been  neglected  the  quantity  4(y  -1)  in  Eq.  (19)  would  be  replaced 

c 

by  2(y  -1)  .  Then  for  Y  =  Y  one  recovers  the  formula  given  in  Zel'dovich 
c  sc 

and  Raizer^. 


In  order  to  include  the  effects  of  thermal  evaporation  we  need  an 
equation  for  the  cavity  temperature  Tc,  which  is  assumed  to  be  uniform 
except  in  the  conduction  front.  The  total  energy  equation  will  be  used  to 
evaluate  the  cavity  pressure  pc  and  by  the  equation  of  state  for  an  ideal 
gas, 

3(M  .  +  M  )  T 

n  *  _ c>ad  evp _ c  (21) 

c  4*[R'3cos2(|e)  +  Rgsin2(|e)]  yc 

In  the  absence  of  evaporative  mass  gain  to  the  cavity,  Mc  acj  is  the  total 
mass  in  the  cavity  under  adiabatic  conditions.  This  mass  is  negligible 


9 


compared  to  that  in  the  swept-up  shell,  but  it  cannot  be  zero  since  the 

cavity  must  have  a  finite  temperature.  The  Taylor-von  Neumann-Sedov 

solution  gives  a  cavity  density  which  depends  only  on  the  self-similar 

variable,  r/R'  .  Hence  the  total  mass  in  the  cavity  increases 
3  s  3 

as  R'  ,  i.e.,  M  ,/p  R"  is  constant,  and  the  porpor tionality  factor 
s  c,ad  o  s 

depends  on  the  artificial  demarcation  in  r/R  separating  the  cavity  from  the 
shell.  We  can  avoid  this  problem  noting  that  for  adiabatic  expansion 


p  3M  ,  Y  -1  P  T 
c  =  .  c,ad  .  s  s  c 

p  =  ^4ttp  R'  T)y  +i  u  x 
S  OSS  c  s 
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from  Eqs.  (11)  and  (21).  During  this  phase,  Pc/ps  is  fixed  by  the  energy 

equation  (Eq.  [20])  and  a  choice  of  the  parameter  Tc/Tg  gives  an  effective 

3 

value  for  M  ,/p  R'  .  We  will  assume  that  this  ratio  remains  constant 

c,ad  os 

after  evaporation  begins.  This  approximation  is  reasonable  in  the  case 
that  9  is  a  small  angle,  resulting  in  a  small  change  to  the  adiabatic 
solution,  and  is  irrelevant  in  the  case  that  the  whole  shell  is  subject  to 
evaporation,  since  MgVp  quickly  exceeds  Mc  a(j. 

Finally,  we  need  to  formulate  a  set  of  equations  for  the  conduction 
front  so  that  Mgvp  and  pevp  can  be  explicitly  evaluated.  We  use  the 
approach  given  in  section  II,  i.e,  the  conduction  front  is  a  thin,  planar 
shell  at  the  inside  edge  of  the  swept-up  shell  and  in  quasi-steady 

A 

equilibrium.  Let  pv  be  the  mass  flow  rate  through  the  conduction  front 
with  v  the  gas  velocity  relative  to  the  front.  Then  setting 


we  have 
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Conservation  of  momentum  through  the  conduction  front  is  simply, 


*2  A  *  2 
pv  +  P  -  Pg  +  Pgvg  , 


(25) 


with  p 


P  (Y  +  1)/(Y  -  1)  since  the  leading  shock  is  taken  everywhere 

OSS 


to  be  a  strong  adiabatic  one.  By  momentum  conservation  the  additional 


pressure  on  the  evaporating  segment  of  the  shell  due  to  the  rocket-like 
effect  is 


evp 


P  v  =  j 
s  s  J 
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This  pressure  prevents  the  conduction  front  from  overrunning  the  leading 
shock. 


For  the  energy  equation  we  balance  the  heat  flux  into  swept-up  shell 
against  the  energy  carried  by  the  mass  flux  off  the  shell.  From  Eq.  (1)  we 
have 
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where  z  =  0  at  the  interface  between  the  conduction  front  and  the  swept-up 
snell  and  increases  toward  the  cavity.  The  generalized  heat  flux  q  in  eq. 
(27)  is  an  approximate  representation  for  the  limitation  on  the  heat  flux 
when  the  electron  mean  free  path  is  on  the  order  of  the  temperature 

length  Up  Under  this  condition  the  electron  distribution  function  is  far 
from  Maxwellian,  the  classical  relation  q  -  -  Kg(T)  $T  breaks  down,  and  the 
maximum,  i.e.,  saturated,  heat  flux  is  that  given  by  the  free-streaming 
limit: 


_,T  .1/2 
1  „  *  a  n  T( — ) 

‘sat  e  m 
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(28) 


The  parameter  oc  has  been  determined  to  be  -  0.06  from  Fokker-Planck 
simulations  and  experiments^  It  is  instructive  to  express  the 
conductivity  coefficient  in  terms  of  the  electron  mean  free  path: 
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and  f(z)  is  a  function  tabulated  by  Braginskii^. 
Eqs.  (28)  and  (29)  gives 


(29) 

Rewriting  Eq.  (27)  via 
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1  “3  .  ,  %  Y  T 
Tpv  +  (pv)  ^r- 


_3/2  ,  e  dTv  ,  r= — 

1  +  f^Z)  Jf 

1  +  “ —  /_£ I— I  'i 
ct/3  (  T  dz 


which  explicitly  displays  the  dependence  of  the  heat  flux  on  the  different 
length  scales. 

C.  Non-dimens ionalization  of  the  Dynamic  Equations 

It  is  convenient  to  write  the  equations  in  dimensionless  form  prior  to 
obtaining  solutions.  For  the  total  mass,  momentum  and  energy  equations  we 
choose  a  fudicial  time  t^  prior  to  which  the  expansion  is  entirely 
adiabatic.  The  value  of  t^  is  determined  by  the  solutions  for  the 
structure  of  the  conduction  front.  At  time  t^  one  has  from  Eqs.  (18), 
(19),  and  (20), 

E  1/5  o/s  9  Ri 

D  s  r  /  ^ \  f.  2/5  n  *  2  1^ 

Ri  5o  fci  ’  Ri  y  i:  ’ 
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We  now  define  the  dimensionless  variables  as 


t  -  Ttr  Rg  =  XRt,  Rg  =  \'*v  pc  =  npc>1 


M  =  a) 
evp 


V1"2'?9’’  Mc.ad  1  W'3  • 


T  s  CTS,  ne  i  vne_s,  pv  ;  jp.R, 


Note  that  the  last  three  definitions  are  written  in  terms  of  the 
evaporating  shell  variables.  They  will  be  used  to  analyze  the  conduction 
front.  Also  recall  from  Eq.  (22)  that  Is  taken  to  be  a  constant 
determined  by  an  initial  choice  of 


<VVi  ■  (VV>1  ■  'c,i 
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where 


.  y  -  1  1  +  Z 
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Equation  (32)  gives  in  turn 


-  M.  cos2(4  9)  X"3,  M  -  M1  sin2(i  9")(X3  -  a>). 

Si  L  Si  L 


^  _  2  1  dX"  2  1  dX 

rs  Yg  +  1  tj  dx  ’  Vs  Yg  +  1  dx 


The  equations  for  energy  (Eq.  [13])  and  momentum  (Eqs.  [14]  and  [15]) 
can  now  be  written,  respectively,  as 
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In  the  last  equation  we  have  used  Eqs  (26)  and  (31).  Equation  (24)  for  the 
rate  of  change  of  the  evaporated  mass  transforms  to 

§?  -  *  *2  It  •  <»> 


Finally  using  Eqs.  (11),  (31),  (32)  the  relation  between  the  cavity 

pressure  and  temperature  (Eq.  [21]  becomes 


25  Ys  -  1  1  +  Zc  ,  SlX'3  +  *  9in  ,,dXN2_ 

n  -  j-  — t ■ ,  i~t~7  t - i - ; — K-rr)  (39) 

Ys  s  X'3cos2(ie)  +  X3sin2(-^6) 


The  sec  of  parameters  is  { y  ,  Y  ,  Z  ,  Z  ,  9 ,  q  }  and  the  set  of 

S  C  S  C  C  f  L 

unknowns  is  { A ,  A',  dA/di,  dA'/dt,  u>,  j}  .  Subject  to  the  initial 
conditions  at  t  =  1  , 


.  dA  dA"  2 

*  '  *•  ar  -  jr  ’  T-  “ 


the  system  of  equations  (36),  (37),  and  (38)  can  be  integrated  forward  in 
time  (f)  ,  using  Eq.  (35)  to  determine  the  cavity  pressure  (n)  ,  once  an 
expression  for  the  evaporated  mass  loss  rate  (j)  has  been  given.  To  do  so 
we  will  need  Eq.  (39)  and  a  solution  of  the  structure  of  the  conduction 
front. 


IV.  RESULTS 

A.  Solution  for  Che  Evaporative  Mass  Loss  Rate 

Let  us  begin  by  multiplying  the  momentum  conservation  Eq.  (25)  by 

2 

p/p  to  obtain 


J2(P*)2 


E.  +  (E— )2  -  (£— )  [i  +  j2  ( 


(p_R  )2 


From  Eq.  (II),  (p  R  )2/p  p  -  (y  -l)/2  ,  and  one  can  now  solve  the 

os  s  s  s 

quadratic  equation  for  the  pressure  ratio: 
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ill  +  j2-^-±  {(1  +  j2-——)2  -  4j2  -3-  2  I  ^i/2j  -(41) 


Since  j  is  limited  to  about  1  and  Y  ■  1.2  from  the  experimental 

o  i  2  8 

analysis'1,  -—J  (y  -  1)  «  1  .  The  positive  sign  is  then  the  proper  choice 

to  satisfy  the  boundary  condition  at  the  front  of  the  conduction  front 

where  u  T/uT  *  1  .  Clearly,  Eq.  (41)  constrains  the  mass  flux  giving 
s  s 

approximately 


(42) 


When  the  equality  holds,  the  maximum  value  for  the  mass  loss  rate  jmx 
corresponds  to  the  maximum  pressure  drop  through  the  conduction  front. 

Note  that  although  the  temperature  at  the  rear  of  the  conduction  front  is 
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that  of  the  cavity  (T  ) ,  the  pressure  at  this  point  (p 2 )  is  larger  than  pc; 

once  the  evaporating  segment  accelerates  relative  to  the  adiabatic 

segment  p  >  p'  giving  p  -  p'  /2  <  p  /2  <  p  .We  surmise  that  the 

s  s  csss 

conduction  front  will  adjust  to  maximize  the  pressure  drop  and  hence  also 

1 2 

j.  In  the  terminology  of  fluid  combustion  ,  the  choice  of  the  positive 
sign  in  Eq.  (41)  means  that  our  conduction  front  is  analogous  to  a 
deflagration  reaction.  The  composite  structure  of  a  conduction  front  and  a 
preceeding  shock  wave  can  be  viewed  as  a  deflagration  initialed  by  a  shock 
i.e.,  a  detonation  wave.  Our  assumption  that  j  31  is  equivalent  to  the 

Chapman- Jouguet  hypothesis. 


The  size  of  the  conduction  front  is  given  by  solving  the  energy 
equation  through  the  conduction  front,  Eq.  (30).  To  do  this  first  scale  z 
by  the  radius  Rg  (z  *  y  Rg)  and  then  combine  Eqs.  (29)  and  (32)  to  get 
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Except  for  the  variation  in  the  Coulomb  logorithm  and  Z,  the  term  in  the 
curly  braces  is  the  ratio  of  the  electron  mean  free  path  in  the  swept-up 


shell  to  R  .  Let  us  then  rewrite  Eq.  (43)  as 
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By  using  this  equation  and  Eqs.  (11)  and  (32),  the  energy  equation  (30)  can 
be  transformed  to 
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and  M  is  the  adiabatic  Mach  number  of  the  flow  relative  to  the  conduction 


front.  We  make  several  approximation  in  order  to  integrate  equation 
(45).  First  we  neglect  the  variation  of  Z  through  the  conduction  front  and 
set  Z  equal  to  Zc  in  Eqs.  (42),  (43),  and  (46).  Second,  we 
set  y  *  Yc  ”  5/3.  It  is  then  reasonable  to  neglect  the 
term  (y  -  l)M2/2  compared  to  unity.  The  Chapman-Jouguet  condition  for  the 
conduction  front  has  M  reaching  a  maximum  of  one  at  the  rear  of  the  front 
and  so  this  term  is  <  1/3.  Finally  we  can  write 


V? 


n  T 
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e,s,  s 


z<1  +  V  p 

zsd  +  Ps 


(47) 


Before  we  noted  that  p/ps  drops  from  1  to  '  1/2  through  the  front  and  from 
Eq.  (41)  most  of  the  change  occurs  near  the  rear.  Thus  we  fix  v?  *  1.  The 
solution  to  Eq.  (45)  is  now  readily  found  to  be 


y  -  (-j^)  iff-  (c5/2  -  i)  -|u3  “  DJ 


(48) 


The  first  property  to  note  from  this  result  is  that  at 
?  ■  Cc»  y  must  be  less  than  unity  in  order  to  fit  the  conduction  front 
within  the  blast  wave  radius  Rg.  Initially  c  has  the  constant 
value  c  .,  while  l  /R  decreases.  The  evaporation  process  can  begin  once 
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By  a  series  of  manipulations  involving  Eqs.  (11)  and  (31)  one  finds 


e,s 


m . 


-1Y3 


Cy3-D 


5  3e  3/5 

O  0 


1  /dX> 


1  2x  e^lnA  Z  2(1  +  Z  j7  (y  +  l)5jp  d/bt_  i ^ /  =>  "Tdr'  ’ 

C  S  S  0  1 


(50) 


We  now  identity  t^  with  the  turn-on  time  for  evaporation. 

At  this  time  A  ■  1,  dX/di  *  2/5  and  by  using  the 

-3  —6 

relation  pQ(gr  cm  )  ■  1.64  x  10  pQ  (Torr)  we  can  convert  Eq.  (50)  to 
practical  units  and  explicitly  define  the  initial  time  as 
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The  corresponding  length  is 
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We  can  further  use  Eq.  (48)  to  determine  a  turn-off  condition  for 

evaporation.  After  evaporation  begins,  the  size  of  the  conduction  front 

yc,  given  by  setting  z,  *  ?c  in  Eq.  (48),  decreases  rapidly.  This  is  due  to 

the  decrease  in  Tc/Tg  as  the  evaporating  shell  segment  accelerates  relative 

to  the  insulated  segment.  Since  there  is  only  a  small  mass  addition  into 

the  cavity  for  small  0,  T  IT'  remains  nearly  constant,  T  /T'  increases,  and 

OS  s  s 

hence  T„/T_  decreases.  Furthermore  l  /R  decreases  as  evidenced  by  the 
c  a  e, s  s 

form  of  Eq.  (50).  At  some  point  y£  can  become  on  the  order  of  the  electron 
mean  free  path  in  the  cavity.  The  latter  quantity  is  given  by 


l  C  1 
e,s  *c 

R  v 

s  c 


We  have  rather  arbitrarily  cut  off  conduction  when  yc  is  twice  this 
value.  The  aspect  represents  the  fact  that  at  this  point  electrons  can 
freely  stream  from  the  cavity  right  into  the  swept-up  shell.  The  resulting 
large  anisotropy  in  the  electron  distribution  is  not  accounted  for  in  the 
present  simple  model  for  saturated  conduction. 


To  summarize,  the  value  of  the  cavity  to  shell  temperature  c  is  given 

c 

by  Eq.  (39).  The  value  of  the  mass  evaporation  rate  (j )  can  then  be 

computed  using  the  equality  in  Eq.  (42).  If  the  size  of  the  conduction 

front  is  less  than  R  (y  <  1),  the  calculated  value  of  j  is  used  in  eq. 

s  c 

(38).  However,  if  the  size  of  the  front  is  less  than  twice  the  electron 

mean  free  path  in  the  cavity,  the  value  of  j  is  set  to  zero.  We  note  from 

-1/2 

eq.  (42)  that  j  a  z,  .  This  is  consistent  with  the  form  of  equation  (45) 

c 

since  the  heat  flux  becomes  more  strongly  saturated  as  the  expansion 

proceeds. 


B.  Solution  for  the  Evolution  of  the  Aneurism 

In  Figure  3  we  present  the  results  for  three  calculations.  For  all 

angels  we  found  that  the  expansion  of  the  insulated  segment  of  the  shell 

2/5 

was  close  to  the  purely  adiabatic  law,  X  =  t  .  This  relation  is  shown  as 

a  solid  line  in  the  graph.  The  position  of  the  evaporating  segment,  i.e., 

the  aneurism,  is  represented  by  the  different  symbols  corresponding  to  the 

different  choices  for  the  parameters  which  are  listed  in  Table  1.  The  time 

and  length  units,  t^  and  Rp  are  calculated  from  Eqs.  (49),  (51),  and  (52) 

for  different  initial  choices  of  5  ,  =*  (T  /T  ) , .  These  units  are 

c,  1  c  s  1 

normalized  to  Eq  *  10  Joules  and  pQ  =  5  Torr,  which  are  typical  values  for 
the  high  pressure  experiments. 

The  calculations  do  not  indicate  a  significant  angle  dependence  of  the 
expansion  of  the  protuberance  for  9  within  a  factor  of  two  of  30°.  However, 
changes  in  the  parameters  Zc  and/or  do  result  in  a  wide  variation  in  the 
shape  of  the  solution  curve  for  Rg/R^ »  as  well  as  changes  in  the  time  and 
length  units. 

V.  SUMMARY  AND  DISCUSSION 

In  this  report  we  have  suggested  a  new  physical  mechanism  for  the 
development  of  the  aneurism  observed  in  the  laser-Hane  experiment.  The 
essential  feature  of  the  model  is  the  formation  of  a  thermal  conduction 

front  on  the  inside  surface  of  the  swept-up  blast  wave  shell.  The 

conduction  front  is  limited  to  a  small  segment  about  the  symmetry  axis  of 
the  incoming  laser  beam  due  to  an  azimuthal  magnetic  field  carried  outward 
by  the  shell.  The  field  is  generated  at  the  target  during  the  laser 
burning.  Within  the  conduction  front  the  thermal  heat  flux  from  the  cavity 
into  the  shell  is  balanced  by  an  evaporative  mass  flow  from  the  shell  back 
into  the  cavity.  The  major  analysis  of  this  report  used  the  thin-shell 
approximation  to  model  the  rocket-like  acceleration  of  the  evaporating 
segment  ahead  of  the  insulated  part  of  the  blast  wave  shell. 

3elow  are  a  number  of  points  on  the  implication  and  limitations  of  the 
present  model  for  the  development  of  the  aneurism. 

I.)  Figure  3  and  Table  I  clearly  shows  that  the  growth  of  the 

protuberance  is  faster  the  smaller  the  value  of  (Tc/Tg)p  This  feature  is 
contrary  to  one's  expectation  in  a  mechanism  where  thermal  conduction  is 


important.  However  this  expectation  is  based  on  the  use  of  the  classical 
thermal  conductivity  coefficient.  Let  us  use  Eq.  (48)  to  determine  j; 


We  can  effectively  keep  the  heat  flux  from  strongly  saturating  by 

maximizing  the  size  of  the  conduction  front,  i.e.,  setting  y  -  1.  Then 

-14/5  C 

as  (l  /R  )  decreases  (roughly  as  t  from  Eq.  [50],  the  conduction 

6  )  S  S 

front  would  become  less  and  less  saturated  and  the  evaporation  rate  would 

5/2 

have  the  expected  behavior  on  the  temperature  ratio  j  *  .  The 

resulting  value  of  j  would  be  very  small  after  t  reaches  a  few.  However, 
this  classical  model  is  inconsistent.  The  conduction  front,  which  would 
extend  over  the  whole  cavity,  has  a  negligible  pressure  drop,  contrary  to 
the  known  relation  p^  -  1/2  pg  .  In  the  pressure  model,  this  pressure  drop 
enforces  the  Chapman-Jouguet  condition  for  the  conduction  front.  This 
condition  is  equivalent  to  the  equality  in  Eq.  (42),  which  gives  rise  to 
the  calculated  dependence  in  T _/T_.  As  a  consequence,  Eq.  (54)  shows  that 
yc  «  1  and  the  heat  flux  is  highly  saturated  through  the  front. 

2. )  As  a  test  of  the  present  model  we  solved  the  system  of  equations 

with  29  *  2it  ,  i.e.,  we  modeled  an  "isothermal"  blast  wave.  Independent  of 
the  parameter  (Tc/Ts ) ^ ,  the  results  were  that  after  the  evaporation  process 
commenced,  the  radius  of  the  shell  Rg  was  1.09  times  larger  than  the  radius 
of  the  Taylor-von  Neumann-Sedov,  purely  adiabatic,  blast  wave,  Ra<j.  This 
compares  favorably  with  the  detailed  self-similar  solution^  for  an 

isothermal  blast  wave  where  Rg  =  1.08  Rad’ 

3. )  We  note  that  the  time  scale  for  turn-on  of  the  evaporation 

process  (Eq.  [51  J)  scales  as  E  ^ ^ •  This  is  consistent  with  the 

experimental  observation  that  the  aneurism  does  not  occur,  or  is  delayed, 
in  high  energy,  low  pressure  shots^.  The  parameters  in  Table  1  for  the 
calculations  in  Figure  3  give  time  scales  for  the  appearance  of  the 
aneurism  that  are  within  a  factor  of  two  of  the  experimental  values. 

However  the  experiment  shows  a  wide  variety  in  the  structure  of  the 
aneurism;  some  grow,  some  grow  then  decay,  while  others  appear  turbulent. 
Undoubtably  there  are  a  large  number  of  factors  which  effect  the  dynamics 


of  the  aneurism.  The  model  proposed  here  has  only  attempted  to  analyze  a 
single  potential  mechanism  in  the  nonlinear  regime. 

4. )  One  of  the  most  serious  drawbacks  of  the  present  calculation  is 
that  it  is  limited  to  one  dimension.  The  two  dimensional  aspects  of 
tangential  flow  in  the  aneurism  and  oblique  shock  could  significantly  alter 
the  evolution  of  the  proturbance  in  our  model. 

5. )  A  second  drawback  in  the  model  is  the  neglect  of  the  evolution  of 

the  topology  and  magnitude  of  the  self-generated  magnetic  field.  We  have 

also  ignored  the  presence  of  an  ambient  magnetic  field.  For  the  parameters 
of  the  present  experiment,  the  ambient  field  would  have  no  dynamic  effect 
on  the  growth  of  the  aneurism;  the  total  magnetic  energy  of  a  800  Gauss 
field  in  a  sphere  of  radius  1  cm  is  *•  0.01  Joules  which  is  far  less  than 

the  several  Joules  Imparted  by  the  laser  to  the  target.  This  conclusion  is 

consistent  with  the  claim  of  the  experimental  group  that  the  aneurism  shows 
no  dramatic  change  whether  the  ambient  field  is  on  or  off  (though  see 
Brecht^).  However,  one  interesting  feature  also  noted  by  the  group  was 
that  the  aneurism  sometimes  appears  to  form  slightly  of  the  symmetry 
axis.  A  possible  explanation  for  this  feature  in  the  context  of  the 
present  model  is  a  follows.  When  the  ambient  field  was  on  it  was  directed 
into  or  out  of  the  plane  depicted  in  Figure  1.  As  the  shell  sweeps  over 
some  of  the  ambient  field,  a  cancellation  with  the  self-generated  field 
could  result  in  a  negligible  magnetic  field  at  an  off-axis  position  thereby 
generating  an  off-axis  aneurism.  This  position  would  depend  on  the 
orientation  of  the  ambient  field.  A  review  of  the  present  data  might  be 
relevant  in  this  matter.  We  also  suggest  that  future  runs  of  the 
experiment  be  made  with  the  ambient  magnetic  field  in  the  plane  of  Figure  1 
to  see  if,  and  how,  the  aneurism  is  affected. 

6. )  A  final  limitation  in  the  model  was  the  approximate  calculation 
for  the  structure  of  the  conduction  front.  We  assumed  a  quasi-steady 
state,  neglecting  variations  in  the  ionization  level  Z,  and  ignored  other 
physical  aspects  such  as  temperature  non-equipartition  between  the  ions  and 
electrons.  Despite  these  limitations  we  can  make  two  general  comments.  If 
evaporation  is  acting  to  drive  the  observed  aneurism  there  would  result  a 
mass  flow  off  the  shell  into  the  region  behind  the  aneurism.  This  mass 
flow  is  much  larger  than  the  flow  off  the  shell  from  the  adiabatic,  i.e., 
not  evaporating,  segments.  One  might  expect  to  find  a  higher  gas  density 


behind  the  aneurism  than  behind  the  other  parts  of  the  blast  wave.  This  is 
in  fact  the  observed  case,  at  least  for  the  electrons  in  the  single  shot 
where  an  inter ferogram  was  obtained^.  This  model  would  also  predict  that 
the  gas  near  the  aneurism  would  have  a  higher  temperature  than  the  gas  in 
the  other  segments  of  the  blast  wave.  This  is  due  to  the  relative 
acceleration  of  the  evaporating  segment  leading  to  a  higher  temperature  in 
the  shell  and  to  the  heating  within  the  conduction  front.  It  might  be 
possible  in  future  experiments  to  measure  the  temperature  near  the  aneurism 
using  spectroscopy. 
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ARLINGTON,  VA  22217 
01 CY  ATTN  CODE  465 
01 CY  ATTN  CODE  461 
01 CY  ATTN  COOE  402 
01 CY  ATTN  COOE  420 
01CY  ATTN  COOE  421 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

01 CY  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND /XPD 
DEPARTMENT  OF  THE  AIR  FORCE 


ENT  AFB 

,  CO 

80912 

01  CY 

ATTN 

XPDQQ 

01  CY 

ATTN 

XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 

HANSCOM 

AFB, 

MA  01731 

01  CY 

ATTN 

OPR 

HAROLD  GARDNER 

01  CY 

ATTN 

LKB 

KENNETH  S. 

W.  CHAMPION 

01  CY 

ATTN 

OPR 

ALVA  T.  STAIR 

01  CY 

ATTN 

PHD 

JURGEN  BUCHAU 

0 1  C  Y 

ATTN 

PHD 

JOHN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  87117 
01 CY  ATTN  SUL 

0 1 C Y  ATTN  CA  ARTHUR  H.  GUENTHER 
01 CY  ATTN  NTYCE  1 LT.  G.  KRAJEI 

AFTAC 

PATRICK  AFB,  FL  32925 
01 CY  ATTN  TF/MAJ  WILEY 
01 CY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WR1GHT-PATTERS0N  AFB,  OH  *5433 
01 CY  ATTN  AAD  WADE  HUNT 
01 CY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20330 
01 CY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 
01 CY  ATTN  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DI V I S ION / YSE A 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01732 
01 CY  ATTN  YSEA 


ELECTRONIC  SYSTEMS  OIVISION/DC 
DEPARTMENT  Of  THE  AIR  FORCE 
HANSCOM  A  F  8 ,  MA  01731 

01 CY  ATTN  DCKC  MAJ  J.C.  CLARK 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  NICO  LIBRARY 

01 CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13441 

01 CY  ATTN  DOC  LIBRARY/TSLO 
0 1 C Y  ATTN  OCSE  V.  COYNE 

SAMSO/SZ 

POST  OFFICE  BOX  92960 
UORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
(SPACE  DEFENSE  SYSTEMS) 

01 CY  ATTN  SZJ 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB,  NB  68113 

01 CY  ATTN  ADWATE  MAJ  BRUCE  BAUER 
01  CY  ATTN  NRT 

01 CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WORLDHAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

01 CY  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AFB,  CA  92409 
(MINUTEMAN) 

01 CY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  MA  01731 

01 CY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

01 CY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

01 CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG8G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

01  CY  ATTN  DOC  CON*FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAURENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

01 CY  ATTN  DOC  CON  FOR  TECH  INFO 
DEPT 

01  CY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
01  CY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 


01  CY 

ATTN 

DOC 

CON 

FOR 

J. 

WOLCOTT 

01  CY 

ATTN 

DOC 

CON 

FOR 

R. 

F.  TASCHEK 

01  CY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

01  C  Y 

ATTN 

DOC 

CON 

FOR 

J  . 

MALIK 

01  CY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

01  CY 

ATTN 

DOC 

CON 

FOR 

J  . 

ZINN 

01  CY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

01  CY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTERVELT 

01  CY 

ATTN 

D. 

SAPPENFIELD 

SANDIA 

LABORATORIES 

P.O.  BOX  5800 

ALBUQUERQUE, 

NM 

87115 

01  CY 

ATTN 

DOC 

CON 

FOR 

W. 

BROWN 

01  CY 

ATTN 

DOC 

CON 

FOR 

A. 

THORNBROUGH 

01  CY 

ATTN 

DOC 

CON 

FOR 

T. 

WRIGHT 

01  CY 

ATTN 

DOC 

CON 

FOR 

D. 

DAHLGREN 

01  CY 

ATTN 

DOC 

CON 

FOR 

3141 

01  CY 

ATTN 

DOC 

CON 

FOR 

SPACE  PROJECT 

DIV 


SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

01  CY  ATTN  DOC  CON  FOR  B.  MURPHEY 
0 1 C Y  ATTN  DOC  CON  FOR  T.  COOK 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  O.C.  20545 

0 1 C Y  ATTN  DOC  CON  DR.  YO  SONG 


OTHER  GOVERNMENT 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  8  INFO 
ADMIN 

BOULDER,  CO  80303 

01 CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01 CY  ATTN  W.  UTLAUT 
0 1 C Y  ATTN  D.  CROMBIE 
01 CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  8  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
01 CY  ATTN  R.  GRUBB 
01 CY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 
P.O.  BOX  92957 


LOS  ANGELES, 

CA 

90009 

01  CY 

ATTN 

I. 

GARFUNKEL 

01  CY 

ATTN 

T. 

SALMI 

01  CY 

ATTN 

V. 

JOSEPHSON 

01  CY 

ATTN 

S. 

BOWER 

01  CY 

ATTN 

D. 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  HA  01803 

01 CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 

01 CY  ATTN  L.  SLOAN 
01 CY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
01 CY  ATTN  J.  WORKMAN 
01 CY  ATTN  C.  PRETTIE 
01 C Y  ATTN  S.  BRECHT 


80EING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

0 1 C Y  ATTN  G.  KEISTER 
01 CY  ATTN  D.  MURRAY 
0 1 C Y  ATTN  G.  HALL 
01 CY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
01CY  ATTN  D.B.  COX 
01 C Y  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
0 1 C Y  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

01 CY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
01 CY  ATTN  H.  LOGSTON 
01 CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wi  IsM  re  Blvd. 

Santa  Monica,  Calif  90401 
01 CY  ATTN  C.B.  GABBARD 
01 CY  ATTN  R.  LELEVIER 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
01  CY  ATTN  J.  ROBERTS 
01  CY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
01 CY  ATTN  M.H.  BORTNER 
SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
01 CY  ATTN  F.  REI3ERT 
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GENERAL  ELECTRIC  TECH  SERVICES 
CO.,  INC. 

HUES 

COURT  STREET 
SYRACUSE,  NY  13201 

0 1 C Y  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01 CY  ATTN  T.N.  OAVIS  (UNCLASS  ONLY) 
01 C Y  ATTN  TECHNICAL  LIBRARY 
01 C Y  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01 C Y  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01 CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01 CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01 C  Y  ATTN  J.M.  AEIN 
01 C Y  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  WOLFARD 
01 C Y  ATTN  JOEL  BENGSTON 

INTL  TEL  8  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01 CY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 

P.O.  BOX  85154 

SAN  DIEGO,  CA  92138 

0 1 C Y  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
0 1  C Y  ATTN  THOMAS  POTEMRA 
0 1  C Y  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01 CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
01 CY  ATTN  DAS  I  AC 
0 1 C Y  ATTN  WARREN  S.  KNAPP 
01 CY  ATTN  WILLIAM  MCNAMARA 
01 CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

01 CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01 CY  ATTN  DEPT  60-12 
01 CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-12 
01 CY  ATTN  W.L.  IMHOF  DEPT  52-12 
01 CY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

01  CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01 CY  ATTN  R.  HEFFNER 


M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  DAVID  M.  TOWLE 
01 CY  ATTN  L.  LOUGHLIN 
0 1 C  Y  ATTN  D.  CLARK 


MCDONNEL  DOUGLAS  CORPORATION 


5301  90LSA  AVENUE 


HUNT 

INGTON  BE 

AC 

H,  CA  92647 

01 

CY 

ATTN 

N. 

HARRIS 

01 

CY 

ATTN 

J  . 

MOULE 

01 

CY 

ATTN 

GEORGE  MROZ 

01 

CY 

ATTN 

W. 

OLSON 

01 

CY 

ATTN 

R. 

W.  HALPRIN 

01 

CY 

ATTN 

TE 

CHNICAL 

LIBRA 

RY 

SERVICES 

MISS 

ION 

RESEARC 

H  CORPORATION 

735 

STA' 

rs  STR 

:eet 

SANTA  BARBARA 

0 

CA  93101 

01 

CY 

ATTN 

P. 

FISCHER 

01 

CY 

ATTN 

W. 

F.  CREVIER 

01 

CY 

ATTN 

ST 

EVEN  L.  GUTS  C 

01 

CY 

ATTN 

R. 

BOGUSCH 

01 

CY 

ATTN 

R. 

HENDRICK 

01 

CY 

ATTN 

RALPH  KILB 

01 

CY 

ATTN 

DAVE  SOWLE 

01 

CY 

ATTN 

F. 

FAJEN 

01 

CY 

ATTN 

M. 

SCHEIBE 

01 

CY 

ATTN 

CONRAD  L.  LONGM 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

0 1 C Y  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01  CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 C  Y  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  8  D  ASSOCIATES 
P.O.  BOX  9695 


01  CY 

ATTN  B.  WHITE 

MARINA 

DEL  REY, 

CA  90291 

01  CY 

ATTN  R.  STAGAT 

01  CY 

ATTN 

FORREST  GILMORE 

01  CY 

ATTN 

WILLIAM  B.  WRIGHT 

MISSION 

RESEARCH  CORP. 

01  CY 

ATTN 

WILLIAM  J.  KARZAS 

1720  RANDOLPH  ROAD,  S.E. 

01  CY 

ATTN 

H. 

ORY 

ALBUQUERQUE,  NEW  MEXICO  37106 

01  CY 

ATTN 

C. 

MACDONALD 

01  CY 

R.  STELLINGWERF 

01  CY 

ATTN 

R. 

TURCO 

01  CY 

M.  ALME 

01  CY 

ATTN 

L. 

DeRANO 

01  CY 

L.  WRIGHT 

01  CY 

ATTN 

W. 

TSAI 

JR. 


MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

0 1 C Y  ATTN  JOHN  MORGANSTSRN 
01 CY  ATTN  G.  HARDING 
01 CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01 CY  ATTN  W.  HALL 
01 C  Y  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  8LV0. 

LOS  ANGELES,  CA  90025 

01 C  Y  ATTN  E.C.  FIELD,  JR. 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
01 CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

0 1  C Y  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01 C Y  ATTN  VINCE  TRAPANI 


* 


SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

0 1 C Y  ATTN  LEWIS  M.  LINSON 
0 1 C  Y  ATTN  DANIEL  A.  HAMLIN 
0 1 C  Y  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01 CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
01  CY  J  .  COCKA'i  - 
0 1 C  Y  E.  HYMAN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 
01 CY  ATTN  J.  CASPER 
01 CY  ATTN  DONALD  NEILSON 
01 CY  ATTN  ALAN  BURNS 
01  CY  ATTN  G.  SMITH 
01 CY  ATTN  R.  TSUNODA 
01 CY  ATTN  DAVID  A.  JOHNSON 
01 C Y  ATTN  WALTER  G.  CHE3NUT 
01 CY  ATTN  CHARLES  L.  RINO 
01 CY  ATTN  WALTER  JAYE 
01 CY  ATTN  J.  VICKREY 
01 CY  ATTN  RAY  L.  LEADA3R AND 
01 CY  ATTN  G.  CARPENTER 
01 CY  ATTN  G.  PRICE 
01 CY  ATTN  R.  LIVINGSTON 
01 CY  ATTN  V.  GONZALES 
01 CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01 CY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01 CY  ATTN  JOHN  ISE,  JR. 

01CY  ATTN  JOEL  GAR9ARIN0 

TRW  DEFENSE  8  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90273 
01 CY  ATTN  R.  K.  PLE3UCH 
01 CY  ATTN  S.  ALTSCHULER 
01 CY  ATTN  D.  DEE 
01 CY  ATTN  0/  STOCKWELL 
SNTF/1575 


VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01 CY  ATTN  W.  RE  ID Y 
01 C Y  ATTN  J.  CARPENTER 
0 1 C  Y  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01 CY  ATTN:  N.  ZABUSKY 


